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ABSTRACT 

The production of o(^) particles in a weakly-coupled theory is believed to 
be non-perturbatively suppressed. I comment on the prospects of (a) establishing 
this rigorously, and (b ) estimating the effect to exponential accuracy semiclassically, 
by discussing two closely-related problems: the large-order behaviour of few-point 
Green functions, and induced excitation in quantum mechanics. Induced tunneling 
in the latter case is exponentially enhanced for frequencies of the order of the barrier 
height. 
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1. Introduction. The issue of baryon- and lepton- number violation in the stan- 
dard electroweak model [1] [2] [3] has brought back to the limelight the limitations 
of perturbation theory. The problem, in a nutshell, is that naive weak-coupling ex- 
pansions cannot be used to estimate processes in which a large number of particles 
is involved. If we denote for instance by N = ^ the number of W bosons produced 
in a high-energy collision, then an expansion of the inclusive cross-section in terms 
of Feynman diagrams corresponds to expanding simultaneously in the gauge cou- 
pling g and in v. This is of no use if one is interested in the region g — > with v 
held fixed and not small. Indeed the leading-order result for this process, whether 
accompanied [2] or not [4] by vacuum tunneling, violates for large v the unitarity 
bounds, and is hence manifestly unreliable. A breakdown of perturbation theory 
of course also occurs in other contexts, such as at high temperature, high densities 
or large external fields. In all these cases it is however obvious that the properties 
of the perturbative vacuum are being modified drastically, so that the theory is 
no more weakly-coupled. In high-energy two-particle collisions on the other hand, 
our difficulty in calculating multiparticle production looks more like a technical 
nuisance, rather than a signal that these processes are unsuppressed. This sounds 
intuitively obvious if one thinks of the time-reversed process: 

how could we send 50 W bosons in an interaction region, and expect only an 
energetic e + e~ pair to emerge? 

There has been much learned debate on this issue, which I could not possibly 
review in this short talk. Here I will focus briefly on two closely-related problems 
where definitive answers are available: the problem of induced excitation in quan- 
tum mechanics, and the large-order behaviour of typical Euclidean Green functions. 
Based on these I will then make a few comments on the prospects (i) of proving 
that multiparticle production in a weakly-coupled theory is non-perturbatively 
small, and (ii) of actually calculating it to exponential accuracy in a semiclassi- 
cal approximation. Given the intuitive evidence for exponential suppression, and 
the ridiculously small low-energy tunneling probability in the electroweak model, 
e ~ 10 , one may wonder whether calculating zero to exponential accuracy 
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is really worth all this effort. I think the answer is yes for two reasons: first, the 
issue is sufficiently important to deserve that we close all loopholes in our intuitive 
arguments. Second, learning how to calculate such non-perturbative phenomena 
may turn out to be academic in the electroweak model, but very interesting in 
other contexts. In this respect simpler models could be more valuable for their 
own sake, rather than as paradigms of multi-VT and -Z production. 



2. Large-order behaviour. We are all well-accustomed to the fact that in field 
theory perturbative expansions are usually divergent. Indeed, a finite radius of 
convergence would be in contradiction with the vacuum instability that typically 
develops for negative square coupling [5]. Consider for definiteness a scalar field in 
d dimensions, 

whose action in appropriate mass units is 

S = Jd d X [±(00)2 + 1^+ ^] , (I) 

and let 

00 



G^Vl Pi,..., PN ) x J29 2n G { n N \pi,..,PN) (2) 

n=0 



be the expansion of a Euclidean iV-point function. For n > 1 we can compute 
the coefficients of the series semiclassically [6] [7]. The idea is to relate them by a 
dispersion integral to the discontinuity of the function on the negative-^ 2 axis, 

o 

G { n N \ Pl ,...,p N ) = ^- J ^discG^Vbi,...,^) (3) 

— oo 

The integral is then evaluated at the saddle point (f) c i(x) = i— f(x), which 
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describes the decay of the vacuum for negative square coupling, with the result 

a N 

G^ipu-.PN) ~„ (2n)^r(n+^)(-ra-^llf( Pj ) . (4) 

3=1 

Here S((p c i) = is the saddle-point action, and we have dropped a momentum- 
conserving 5-function. Note that a weak-coupling expansion of the discontinuity is 
justified in the large-n limit, provided N and pj are all kept finite and fixed. 

The rapid factorial growth of the coefficients, typical of bosonic field theory, 
makes the series (2) diverge for any g. A closer look at Feynman diagrams in fact 
reveals that these large contributions at high orders arise when n ~ o(^) vertices 
are concentrated in a real-space region of size ~ o(l) *. From the functional-integral 
point of view these contributions probe the region of large fields, or of large virtual 
non-linear waves. It is because such waves are not described adequately, when one 
expands around the usual vacuum, that perturbation theory breaks down. Intu- 
itively we of course expect that large fluctuations, though perhaps hard to calculate 
precisely, are nevertheless exponentially suppressed and thus often negligible. This 
expectation is in practice confirmed by the enormous success of QED: keeping for 
instance three terms in the expansion of the anomalous magnetic moment of the 
electron, we find agreement with experiment to better than seven significant digits 
[9]. The more precise mathematical statement one would like to make is that the 
perturbative series is asymptotic to some rigorously defined function, and that the 
semiclassical estimate, eq.(4), can be converted into a uniform bound for the series 
remainder: 

Rem n (G) = \G^ N \g 2 ) - £<?, (J V| < n! g 2n a n , (5) 

where lim n _ >00 (a n ) _1 / n = a is a finite constant, and g 2 lies in some finite interval 
(0, g 2 ) on the positive real axis. From the above bound and from Stirling's formula, 



* Fermionic theories fare in this respect better thanks to the Pauli exclusion principle. In 
closed-string theory, on the other hand, the divergence of perturbative expansions appears 
to be even more severe [8]. 
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n\ ~ \phx[n + l) n+ 2e~ n_1 , we could conclude that in the g — > limit 

min Rem n (G) ^ exp^ — ^ j > (6) 

so that the divergent high-order contributions indeed sum up into a controllable 
non-perturbatively small ambiguity . 

Proving this statement has been one of the aims of formal (constructive) field 
theory [10]. For the Euclidean scalar theory, eq.(l), in the superrenormalizable 
(d < 4) domain, the stronger result of Borel summability [11] ensures in particular 
the existence of the above bound K Furthermore, the large-order semiclassical 
estimate can be shown to correctly predict the relevant singularity of the Borel 
transform [13], so that the bound is optimized for a = a. The extension of these 
results to (%) the Minkowski region in field theory, (%%) the double- well potential in 
quantum mechanics, or (Hi) the d — 4 renormalizable theory is not straightforward. 
Controlling non-perturbative effects in the first case is hard because bounds do not 
continue analytically. Borel summability has nevertheless been established for the 
on-shell four-point function, but only below the particle-production threshold [14]. 
In the case of the double- well potential of quantum mechanics, obtained by flipping 
the sign of the quadratic term in (1), there is a hindrance to Borel summability due 
to the presence of real instantons. Indeed the semiclassical large-order analysis [15] 
predicts a singularity of the Borel transform at the same distance, a = |, from the 
origin as in the single-well case, but lying on the positive real axis. This singularity 
corresponds to the action of an infinitely-separated instanton-antiinstanton pair, 
and is expected to control the non-perturbative ambiguities of the asymptotic 
expansions in the trivial and one-instanton sectors [7], but rigorous bounds of the 
type (5) have not to my knowledge been established [16]. In the third case of 
the 0| model, Borel summability is obstracted by the so-called renormalons [17], 

* The simple statement of asymptoticity of the series says nothing about the size of this 
ambiguity, which could for instance have been as large as exp(— ). 

J Borel summability in fact allows a reconstruction of the function from its series with arbi- 
trary precision, by using appropriate conformal mappings [12] [7] . 
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which are a warning that non-perturbative ambiguities could a priori be affected by 
renormalization. Since, however, the theory most probably does not exist, there 
is little point in worrying about bounds (5). Needless to say, finally, that non- 
perturbative effects in four-dimensional gauge theories are still beyond rigorous 
technical control [10], since in addition to renormalization and the presence of real 
instantons one must also face the problems of gauge fixing and scale invariance. 

There are three lessons to retain from this blitz review of large-order behaviour: 
first that naive perturbation theory around the vacuum breaks down when one tries 
to estimate the contribution of large fields, or of (real or virtual) non-linear waves. 
Second that such contributions to few-particle processes are expected to stay ex- 
ponentially small, and can be estimated to exponential accuracy semi-classically. 
And third that a real proof of exponential suppression requires a non-perturbative 
control of the vacuum and takes us into the rough territory of constructive field 
theory. These comments should be kept in mind when moving on to the problem 
of multi-particle production, which in a vague sense is the square root of the large- 
order problem: indeed, large-field fluctuations are in this case created from, but 
do not have to disappear back to a few-particle state. 



3. Multi-leg functions in Quantum Mechanics. The quantum mechanical ana- 
log of multiparticle production is the induced excitation of an anharmonic oscillator 
under the action of a weak but very energetic external force [18] [19] [20]. We con- 
centrate first on the single- well potential with action given by eq.(l). The ampli- 
tude of interest is proportional to the matrix element (O|0|iV), where N = ^> 1 
is the level of the excited final state and E = is its energy. To motivate this 
scaling of parameters consider the weak-coupling expansion of the energy, 

1 3o 2 

E = (N + -) + -f-(2N 2 + 2N + 1) 

2 16 4 (7«) 

- — (34iV 3 + 51iV 2 + 59N + 21) + ... 
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which can be clearly reorganized as follows: 

6 = e(v)+g 2 e 1 (is)+g 4 e 2 (is) + ... (7b) 

The effective expansion parameter in (7b) is Planck's constant, as would be made 
obvious if we were to restore all units in our equations. In particular, the function 
e(v) = v + |z/ 2 + y|g^ 3 + ... is the inverse of the integrated classical density of 
states, 

*> = 1 7 dx )/'('- t-t) (8) 

o 

where (f) = x/g defines the rescaled position variable, and x(e) = \fy/\ + 4e — 1 
is the rescaled classical turning point. This relation between energy and level 
illustrates that, for coherent states of many quanta, it is the classical but not the 
naive weak-coupling approximation that is adequate. 

Let us then consider an analogous semiclassical expansion of the matrix- 
element (N'\x M \N) , when all quantum numbers N = jp, N' = ^, and M = ^ 

are large \ 

(N'\x M \N) =exp( ^F(z/,z/,/i) + Fi(z/,z/, /U ) + ...) . (9) 

The form of this expansion is based on our expectation that in some range of 
parameters the overlap integral should be non-perturbatively suppressed. We are 
of course ultimately interested in the matrix element (O|0|iV), since we want to 
mimic the effect of a few-particle initial state. We must therefore hope that this 
matrix element can be obtained, at least to exponential accuracy, by taking the 

f Because of reflection symmetry N + N' + M must be even. 
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/j,, v' — > limit, i.e. in formulae: 

lim £ 2 log(O|0|iV) = F{y) , 10) 

where 

F(v)= lip F{v,i/,n) . (11) 
The simple unitarity relation 

oo 

(O|0 2 |O) = ^ |(O|0|m)| 2 (12) 

m=0 

implies in particular that |(O|0|iV)| 2 < (O|0 2 |O) ~ \ + o(g 2 ) , so that the func- 
tion F(u) must clearly stay non-positive. The first term of a naive weak-coupling 
expansion violates this unitarity bound, because of the same factorial growth of 
graphs which is responsible for the large-order divergence. A precise calculation in 
fact gives [21]* 



mm = 




(13a) 



from which we find 

F(v) = ~(l + log-) + o(v 2 logv) . (136) 

_ IS 

The Born approximation thus hits the unitarity bound at v — 16e, at which point it 
is manifestly unreliable. In reality, however, the matrix element stays exponentially 
small for all finite v. This can be established explicitly both (a) by a semiclassical 
calculation [22] [18] [23], and (b) by deriving rigorous upper bounds [20], as I will 
now briefly explain: 

* The calculation is identical to the Born approximation for the production of N scalar 
particles at threshold, up to a normalization -^j, and an extra factor of \/2 per external 

leg appearing in the LSZ reduction of quantum mechanics [20]. 
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(a) The semi-classical calculation of the overlap integral J x M ^ N ,^ n for N ^> 
N' ^> 1 was suggested a long time ago by Landau [22]. The idea is to deform the 
integration contour in the complex-a; plane away from the classical turning points, 
use a WKB approximation for the wavefunctions \I/jv an d ^N', an d evaluate the 
resulting integral at its saddle point. The full details of a carefull calculation have 
only recently been completed [23], and confirm Landau's result [18, 22] in the 
H — > limit: 



F(e,e')=- J 




dx 

T-U) 

(14) 



where K is the complete elliptic integral. The result is here given as a function of 
energies, but can be also expressed in terms of the levels v and u' via the classical 
relation (8). As was already anticipated, F has a smooth v' — > limit, which can be 
used to estimate (0\(f)\N) to exponential accuracy. Furthermore F(u) is a monotone 
decreasing function of energy or level, so that the matrix element for finite v is 
indeed exponentially small. Notice also that by using the asymptotic expansion 
K(l — 5) ~ /o(?(-y / /|) + o(5 log5), one can recover the Born result, eq.(136). 

(b) The more rigorous proof of exponential suppression [20] makes use of an 
exact recursion relation between matrix elements of powers of 0. Defining the 
recursion coefficients (N'\<f) M \N} = n { ^ N ' ] (N'\<f> M+2 \N) one finds' 1 ": 

2(M+l\ 

j£ 9 { 2 ) f 15) 

M Q{ M A )n M -2nM-A + 2{ M 2 ){E + E')n M ^ + {{E-E'Y-M^' 



f This equation allows a fast numerical evaluation of the low-lying energy levels and wave- 
functions with very high precision [24] . Notice also that in the semiclassical limit it reduces 
to 

^rf + M 2 (e + e') ^ + [(e e' f - M 2 ] r, - £ = 

where — \logr\ = f^-(e, e', m)- It can therefore be used to determine the /^-dependence of the 
function F completely. 
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where (?) are the binomial coefficients, and E and E' are the energies at the N and 
N' levels. This is the analog of Schwinger- Dyson equations, which in field theory 
are instrumental for proving the uniform remainder bounds (5). The problem we 
are here facing is very similar: indeed, since the first N terms in the asymptotic 
power series for (O|0|iV) are zero, we have for all n < N 

Rem n ((O|0|iV>) = (O|0|iV>, (16) 

so that bounding the series remainder up to this order would automatically put 
bounds on the matrix element itself. Skipping further details on the derivation of 
these bounds, let me simply point out that they finally take the form 

F(e) < min B(e,n , (17) 
~~ o<c<f 

where £/g 2 plays the role of the order of the expansion. This should be chosen 
appropriately so as to optimize the bound B, making sure in particular that the 
coupling-constant suppression is not overwelmed by factorial growth. The optimal 
bound has been shown to decrease monotonically with energy [20], like the semi- 
classical estimate eq.(14). Proving that the two actually coincide would require 
some more work. 

The double- well potential, V(<j>) == fg (<j) 2 - £)' \ can be treated with similar 
techniques. Induced excitation, whether accompanied or not by quantum tunnel- 
ing, can be again bounded by an exponential envelope whose decay with energy is 
monotonic [20]. The leading-order instanton calculation [19], on the other hand, 
predicts that induced tunneling grows fast with energy in the low-energy region. 
To be more precise, the transition amplitude from the ground state in one well 
to the iVth excited state in the other is given to exponential accuracy, and in the 
limit of small is, by 

4 v 16 
F(u) = -- + -(1 + log— ) + subleading . (18) 

This rapid initial growth with energy is the result of two competing effects: the 
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enhancement of spontaneous tunneling, starting from some optimal jumping state, 
is partially compensated by the fast-falling probability of exciting this latter from 
the vacuum [25] . The presence of extra real turning points complicates in this case 
the WKB analysis, and the calculation of the full function F{y) has not yet been 
completed [23]. The exponential envelope tells us, however, that it has to reach a 
maximum at some finite distance below the zero-axis, before dropping indefinitely 
in the v — > oo limit. This means that induced tunneling has a resonance of 
exponential proportions, which could, for instance, be observable in semi-conductor 
physics [26]. 

4. Comments on Field Theory. I will conclude with a few telegraphic comments 
on the prospects of extending the results reviewed in the previous two sections to 
multi-particle production in field theory. We saw that the large-order behaviour 
shows no qualitative difference as one passes from (d — 1) quantum mechanics, to 
(d = 2, 3) superrenormalizable Euclidean theories. It should therefore be possible 
to prove that, like induced excitation in the former, multiparticle production in the 
latter is exponentially suppressed. Repeated use of the Schwinger-Dyson equations 
could indeed lead to bounds for Euclidean multi-leg Green functions, but there is 
at present no obvious strategy on how to extend such bounds to the Minkowski 
region. Closely related are efforts to establish bounds by exploiting the relation 
between the forward elastic amplitude and the total inclusive cross-section [27] 

o-indiVs) = 7= 1 = knAei{y/s) . (19) 
V( s - 4 ) s 

This is the analog of the quantum-mechanical unitarity relation, eq.(12). Taking 
the remainder at order N, we can express the inclusive cross section for producing 
at least iV particles in the following form: 

o- n >N(Vs) = 7= Im Rem N (A e i) + Rem N (<J n <jv) • (20) 

V(s-4)s 

Rigorous bounds on the large-order behaviour of the right-hand side could thus 
be used to establish exponential suppression of the left-hand side. Though very 
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plausible, this strategy stumbles again on the fact that little can be said rigorously 
about the large-order behaviour of perturbation theory deep in the Minkowskian 
region. 

In view of these difficulties, attempting to calculate the process semi- 
classically^ remains, I believe, the most promising prospect. This prospect has 
received a boost recently thanks to a clever suggestion by Rubakov and Tinyakov 
[29], who proposed to first calculate a quantity which appears to have a manifest 
semi-classical expansion. This quantity is the inclusive cross-section for an ensem- 
ble of N' = incoming particles, distributed randomly in phase space, but with 
fixed total energy in the center-of-mass frame. Assuming we can calculate it, we 
may then hope to take the v' — > limit, so as to recover the leading exponen- 
tial estimate for processes with only few particles in the initial state. This has 
been illustrated explicitly in our discussion of the single-well potential in quantum 
mechanics. There is, furthermore, some preliminary evidence for the smoothness 
of this limit in the one-instanton sector of the standard electroweak model [30]. 
Whether a useful solution to this complicated saddle-point problem can, however, 
be found remains to be seen. 

I have benefited from discussions with K. Gawedzki, G. Grunberg, J. Lascoux, 
A. Mueller, E. Mottola, P. Tinyakov E. Papantonopoulos and R. Seneor. I am 
particularly indebted to V. Rivasseau and J. Magnen for many patient explanations 
of constructive field theory methods. 



f That a part of this problem, i.e. the calculation of final-state corrections, admits a semi- 
classical expansion has been known for a while [28] [3] . 
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